Spin current rectification from a spin-biased quantum dot 
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We compute the spin current rectification coefficient of a non-equilibrium quantum dot subject to 
a spin bias and an ac charge bias with small amplitude. As a function of the position of the resonant 
level the spin current rectification coefficient shows a set of three peaks around the gate voltage at 
which the resonant or the upper level of the dot is in the vicinity of the equilibrium Fermi level in 
the leads. The peak heights can be related to the average number of the quantum dot electrons. We 
discuss the frequency dependence of the spin current rectification coefficient as well and emphasize 
the effects of the photon-assisted spin transport through the dot. 



Recent advancement in spintronics 043] and quantum 
computing Q results from an increasing interest in 
studying the spin-polarized transport through nanostruc- 
tures [5|. It is advantageous to manipulate spins in- 
stead of charges in materials because spins have longer 
coherence length as well as relaxation time since, gen- 
erally, the spin is more weakly coupled to the environ- 
ment then the charge. As a part of the research effort in 
this field, the spin-polarized transport through a quan- 
tum dot is extensively investigated. It has been shown 
that a quantum dot with Zeeman splitted states oper- 
ates as a phase-coherent spin pump jjj. Spin-polarized 
current has been demonstrated in a quantum point con- 
tact @ and a Coulomb-blockaded quantum dot @. The 
problem that stands at the forefront in this context is 
to devise means of purely electrical control of the spin- 
polarized current. Recently, it has been proposed to cou- 
ple the quantum dot to leads which have spin-dependent 
electrochemical potentials, the so-called spin bias [Io| . 
The spin bias can be applied instead of or in addition 
to a charge bias to quantum dots. It is experimentally 
realized as a spin accumulation at biased contacts be- 
tween ferromagnetic and nonmagnetic materials or when 
a semiconductor lead is illuminated by a circularly polar- 
ized light fllj . The dc transport properties of spin-biased 
quantum dots have been investigated over the last sev- 
eral years 12]. The transport through a quantum dot 
subject to both constant spin bias and time-dependent 
charge bias has recently been addressed as well [13, [l4[ . 

Future applications of spintronic devices will require 
understanding of the time-dependent spin (and charge) 
transport through them. The work on studying the time- 
dependent charge transport through nanoscopic devices 
(quantum dots in particular) has been an active field of 
research over the last 20 years [l5| . Similar type of inves- 
tigation have to be performed concerning the spin trans- 
port in presence of time-dependent fields as well. 

In a previous work 1J] we studied the linear response 
of a quantum dot, subject to a spin bias, to an external 
time-dependent charge ac bias. We examined the be- 
haviour of the linear response spin admittance as a func- 
tion of the gate voltage and the frequency of the applied 
ac field. The linear response admittance is calculated 
by expanding the total current through the dot to linear 
term in the amplitude of the ac field. In this work we ex- 



tend the investigation to take into account the non-linear 
terms and calculate the spin rectification coefficient of the 
dot. When the current- voltage characteristic of a system 
is non-linear or is not odd with respect to the applied 
voltage a dc current can flow through the system when 
an ac bias acts upon it (current rectification). This effect 
has various applications in electronic devices. 

Spin current rectification has been studied, for exam- 
ple, in the case of molecular fl6| and quantum wires 
17]. Spin-sensitive rectified current has been experimen- 
tally observed in scanning tunneling microscope [l8| . Re- 
cently, spin rectification has been studied in nonitinerant 
one-dimensional quantum spin chains [l!|. In this work 
we consider a quantum dot coupled to two leads with 
spin-dependent electrochemical potentials (spin bias) and 
additional time-dependent charge field. We assume that 
the amplitude of the ac field as small and develop the 
tunneling current to second order with respect to the am- 
plitude of the ac field. The spin and charge rectification 
coefficients are derived from this expansion. 

The Hamiltonian of the dot, coupled to two non- 
magnetic leads, which are considered as ideal reservoirs, 
with external time-dependent charge bias applied to the 
dot is given by 



H = Hn 



Hp 



(1) 



Hq is the Hamiltonian of the system, consisting of a quan- 
tum dot, coupled to the leads 



H = e c ^ c t c <y + Un^n^ + ^ £pkap klJ a [jko 

a ka/3 

+ y^^TpkC+apka + H.c), 



(2) 



where c a {c^) is the annihilation (creation) operator for 
an electron with spin a =t, 4, (+, — ) in the dot, e c is the 
energy level and U is the Coulomb repulsion energy be- 
tween two electrons with opposite spins, simultaneously 
residing in the dot. apka(o.p ka ) is the annihilation (cre- 
ation) operator for an electron in the lead j3 — L,R 
with quasimomentum k and spin a. We assume that 
the electrochemical potentials in the leads are spin de- 
pendent, that is, a spin bias is applied to the system: 
ItLa =£f + {eV c + o-eV s )/2, [i Ra = e F - {eV c + <reV s )/2, 
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where V c is the charge voltage, V s is the spin voltage, and 
Sp is the equilibrium Fermi level energy. In writing these 
expressions we assume that the spin is quantized along 
the same axis in both leads. Finally, Tpk is the tunnel- 
ing matrix element. The Hamiltonian H ext describes the 
external time-dependent perturbation. We consider the 
external charge ac bias to be applied to the quantum dot 
through the (central) gate electrode and the correspond- 
ing Hamiltonian takes the form 



G 

H ext = -e-^-u{t)Y,4{ty.{t). 



(3) 



arate equilibrium electron distribution. The leads re- 
main in steady state in the presence of the external ac 
signal because it is applied only to the quantum dot. 
Correspondingly, the leads' GF are given by the ex- 
pressions for nonintcracting equilibrium electron systems 
[221 ] . The steady-state distribution GF for the dot elec- 
trons does not depend on the " center-of-mass" time T 
and is given by G a< (cu) — G ar (uj)T,o a - < (uj)G CTa (uj) where 
£ <t<M = Efe/3 \ T /3k\ 2 A l3< (k,uj) and A f;< (k,uj) is the 
distribution GF for the lead (3. The average number of 
quantum dot electrons < n a > is to be obtained by self- 
consistently solving the equations 



In a model of the dot with geometric capacitancies [20j , 
CG)Ci>CR are the capacitancies of the gate, left, and 
right electrodes, respectively. Cs = Cq + Cl + Cr is the 
total capacitance of the dot. The ac charge bias is u(t) = 
uq exp (ifit) + c.c. Note that throughout the remainder 
of this text (unless explicitly written) we have omitted 
Planck's constant fi. 

Next, we introduce the Keldysh Green's functions 
(GF) of the dot electrons (similar definitions hold for 
the lead electrons as well). There are three types of GF 
in the Keldysh formulation: the retarded (advanced) and 
the distribution functions [2l|. The retarded GF is de- 
fined as G ar (t) = -i6(t) ({c a {t),c+(Q)}) and the distri- 
bution function is G a< (ti,t2) — \ {c^ (i,2)c a (ti)) . The 
curly brackets denote anticommutator and the averaging 
is with respect to the full time-dependent Hamiltonian of 
the dot. In the Keldysh technique the distribution GF 
depends on two times - the "relative" time t = t\ — t% 
and the "center-of-mass" time T = (t\ +t-i)/2, the re- 
tarded (advanced) Green's functions depend only on a 
single time variable. 

We approximate the retarded GF of the dot electrons 
with the simplest expression which captures the basic 
physics and is applicable for temperatures higher than 
the characteristic Kondo temperature. The GF is ob- 
tained by applying a mean-field approximation in the 
equation of motion for it and the result is 
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(4) 



Here, < n a > is the average number of electrons with 
spin a in the dot and 27 is the total elastic level width, 
7 = 1l + in and 79 = tt J2k \Tpk\ 2 S{u - sp k ). In the 
usually employed approximation of energy-independent 
tunneling matrix elements and with a broad flat density 
of states for the lead electrons 7^ are constants. Quali- 
tatively, the spectrum of the quantum dot electrons with 
spin a consists of a lower (resonant) level with energy e c 
and relative weight 1— < n_ CT > and an upper level with 
energy e c + U with relative weight < n_ CT >. 

We calculate the distribution GF supposing that the 
relaxation processes in the leads are much faster as com- 
pared to those in the dot. With a spin bias applied 
to the quantum dot the leads are not in equilibrium 
but we can assume that each spin subband has a sep- 



< n a >- 



da; 
2^ 



Im G a< {io). 



(5) 



The total current Ip°*,/3 = L,R which flows through 
the left/right barrier of the dot is given in the model of 
geometric capacitancies by [20|, |23| 



rtot _ j 
1 L/Ra - 1 L/R 



c 



R/L + Cq Ch/R t 
~ lL/Ra ~ 1 



Cv 



Gy 



l R/Lo 



(6) 



where Ilir are the displacement currents through the 
parasitic capacitancies. They do not depend on the spin 
bias and are linear with respect to the applied ac field 
and we shall not give the explicit expressions for them. 
The tunneling current through the barrier /3 = L, R is 
given by 




[T pkC +(t) apka (t)~H.c.]S(t)) (7) 



with S(t) = f p cxp[-i fdt'H ext (t')]. f p is the closed- 

v 

path time ordering operator [2l| and the brackets (} de- 
note averaging with respect to the time independent part 
of the dot Hamiltonian (namely, Hq). When the dot is 
subject to charge and/or spin bias both the charge I c and 
the spin I s currents can flow. They are defined as follows 
Ic = h + h, Is = {If - h)/e where l a = J2p Ip<y- 

The nonlinear response coefficient is calculated by ex- 
panding the tunneling current in powers of the external 
ac voltage bias u(t). The term quadratic in the amplitude 
Uq is obtained as 



*S(t) 



ia 2 e 3 
2h 3 



dt'dt v u(t')u(t v ) x 



T L kcrC+(t)aLkcr{t) - H.C 



4' 7 (*')<V 7 V)fh E c t'Y (*" ) c -v (*")%' ) (8) 

cr'7 / 

where the subscripts 7, 7' = ± refer to the positive (+) 
or negative (-) branch of the closed-time path, £,+ —£,- — 
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?7+ = — T]- = 1 'Til . The time- variable of a physical quan- 
tity (the current) is always on the positive time branch. 
The coefficient a = Cg/C^. The reader is referred to 
[24j for an introduction of the way similar expressions are 
dealt with. Next we perform Fourier transformation and 
obtain the tunneling rectification coefficient of the dot. 
ft is defined as the ratio of the zero-frequency Fourier 

(2) 

transform of I L J. to the square of the amplitude of the 

external time-dependent field = 7^(0) /mq. Similar 
definition in which the Fourier transform at frequency 2S1 
is involved gives the second-harmonic generation coeffi- 
cient. 

The explicit expression for the tunneling rectification 
coefficient can be cast into the form 



a Lc , = -47ri-^3-7 L J duj 

— oc 

{}Lo{u)[G aa (Lu)G aa {u + Cl)G aa (uj) - c.c] - 
[G ar {Lo)G a< {uj + n)G aa {Lo) 
+G a< {u)[G <jr {oj)G ar {uj + Cl) + c.c.}} + 

(n-»-n)>. (9) 

Here, is the Fermi-Dirac distribution function for 

the electrons in the left lead with the appropriate spin- 
dependent chemical potential. The notation (Cl —Cl) 
means that one must add the same terms as those explic- 
itly written with the substitution of Cl with —Cl. Note 
that this expression is only valid in the approximation 
used to write down the retarded GF for the dot elec- 
trons. The total rectification coefficient a L ot a is calculated 
from Eq. ©. 

The numerical results are obtained with a set of values 
for the model parameters that are appropriate for the real 
systems: jl — 1r = 7, U = 2O7, temperature is taken 
to be O.O57. Throughout the remainder of the text we 
concentrate on the case of pure spin bias applied to the 
dot - V c = 0. Let us only mention that in this case the 
charge current rectification coefficient (measured in units 
of a 2 e 3 /h 3 ) is generally two orders of magnitude smaller 
than the spin one (in units of a 2 e 2 /h 3 ). Consequently, 
the rectified charge current through the dot is very small. 

In Fig. I we present the results for the dependence of 
the total spin rectification coefficient aL S — i a L°\ — a L°l)l e 
on the position of the quantum dot resonant level, that is, 
on the gate voltage applied to the dot. The behaviour of 
Qls for a dot with applied spin bias is qualitatively sim- 
ilar to the behaviour of the charge current rectification 
coefficient for a dot with applied charge dc bias [2CJ, |25j . 
There are two characteristic three-peak structures (one 
central and two side peaks) for positions of the resonant 
level energy close to and —U. They correspond to the 
cases when the resonant or the upper level is near Ep, cor- 
respondingly, the spin is transferred through the dot via 
the resonant or via the upper level. The two side peaks 
are of non-equal heights. The upper side peak (the one 
at higher gate voltage) for the upper level and the lower 



side peak for the resonant level are of suppressed height. 
This behaviour is to be attributed to the dependence of 
the relative weights of both peaks in the density of states 
on the average number of dot electrons < n a > [20| . 
Consider the two side peaks corresponding to the reso- 
nant (higher) level. Each GF in the expression for the 
tunneling rectification coefficient contributes a factor 
f- < 71/ — fj > (< tl — fj >), hence a L<J ~ (1— < n_ CT >) 3 
in the vicinity of e c ~ Ep — 0. Similarly, in the vicinity 
of the upper level aLa ^< ri-a > 3 - The average num- 
ber of dot electrons < n a > diminishes with increasing 
e c (there are more empty states in the leads available 
for the tunneling electron). The spin current is given as 
a difference between the up-spin electron current (which 
depends on < n± >) and the down-spin electron current 
(which depends on < n<|- >) and when a charge bias is 
not applied < >=< >. Hence, in this case the 
heights of the side peaks in the spin current can directly 
be related to the average number of dot electrons as in 
the case of pure charge current [2(| • 

In the case of charge current rectification the dis- 
tance between the two side peaks is 2CI [2(| In 
the present case, however, the photon-assisted tunnel- 
ing involves transitions from the dot states (for a given 
spin) to the spin-dependent electrochemical potentials 
Pl/Rct which are separated by the spin bias [recall that 
eV s = ^Lt ~ A 1 14 = ~(PRt ~ MRJ-)]- Tlie peaks com- 
ing from the up-spin and from the down-spin electrons 
are displaced by =peV s /2. Consequently, the distance be- 
tween the peaks is ~ 20 + eV s . Actually, the side peaks 
are somewhat more distantly separated. As has been 
pointed out previously [2(| (in the case of charge trans- 
port through a quantum dot) the effect of the thermal 
fluctuations is to shift the positions of the peaks to higher 
energies /frequencies. 

Let us now discuss the frequency behaviour of the spin 
current rectification coefficient. The results are presented 
on Figs. 2 and 3. The presented data are distinguished 
by the positions of the resonant level relative to the spin- 
dependent electrochemical potentials of the leads. When 
both resonant and the upper level are below Hl/r 17 (Fig- 
2, lower panel) the spin current rectification coefficient is 
generally a decreasing function of the frequency because 
with increasing CI the electrons cannot follow the applied 
field and the current diminishes. However, there are two 
resonant-like features in the vicinity of CI = \ep — e c \ and 
Cl = \ep — e c — U\. These features are clear evidence 
for a photon-assisted spin transport through the dot - 
the electrons absorb photons in order to tunnel through 
the dot and, as a result, a spin is transferred between 
the leads as a spin current. Similar behaviour is also 
observed in the case of the the so-called empty orbital 
regime in equilibrium when the resonant level position is 
above all electrochemical potentials of the leads and there 
is resonant-like enhancement of the transport coefficient 
for frequencies Ct — e c and Cl = e c + U . 

The upper panel of figure 2 presents the results for 
aL S when the resonant level is below and the upper level 
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is above Ef- In this case ol s changes sign at Q ~ 57 
and O ~ I57. Let us assume that initially the resonant 
level is occupied by a down-spin electron. At frequency 
f2 = \e c \+eV s /2 ~ 57 opens channel for inelastic transfer 
of down-spin electrons through the dot - a down-spin 
electron can absorb a photon with appropriate energy 
in order to tunnel from the resonant level to the right 
lead. The transfer of the up-spin electrons is Coulomb- 
blockaded for these frequencies. For higher frequency 
f2 = e c + U — eV s /2 the up-spin electrons can inelastically 
tunnel through the dot by absorbing a photon, hence, 
the sign of the spin current changes. Analogously, if the 
dot is initially occupied by an up-spin electron at lower 
frequencies the up-spin electrons can inelastically tunnel 
from the dot to the left lead and at higher frequencies 
the down-spin electrons can tunnel from the dot to the 
right lead giving the same sign of the spin current as in 
the previous case. 

In Fig. 3 we show the behavior of a£ S (fl) in the case 
when the resonant (or the upper) level is placed between 
the electrochemical potentials of the leads fi^ < e c < 
fiLf (or fiLi < e c + U < As evidenced the spin cur- 

rent is negative and there is resonant-like enhancement 
for 51 ~ 2O7. Consider the hrst case. If the resonant level 
is occupied by an up-spin electron then, at frequency 
il ~ e c + U + eV s /2 a down-spin electron inelastically 
tunnels out of the dot to the right lead, hence, the down- 
spin current flows from the left to the right lead. If the 
resonant level is occupied by a down-spin electron the up- 
spin current flows from the right to the left lead. Similar 
reasoning applies in the case when /xlj, < e c + U < /iLt- 
Apparently, if the sign of the spin bias is reversed the 
sign of the rectified spin current will also be reversed. 

To the best of our knowledge, there are no experimen- 



tal investigations of the quantum dot setup that we stud- 
ied in this work. But the approximations we have made 
are the same as those usually employed to discuss the 
spin and/or charge transport through a quantum dot. 
We have considered the lead electrons as noninteracting 
particles. This means that we do not self-consistently de- 
termine the tunneling matrix elements and single-particle 
energies but simply take the corresponding quantities as 
input parameters. Also, we do not take into account the 
effect of the time-dependent potential on the leads elec- 
trons, that is, we assume that the ac bias is completely 
screened in the leads. Evidently, the frequency of the ac 
field must be smaller than the plasma frequency of the 
leads. Typically, the plasma frequency of the leads is 
of the order of tens of GHz. Hence, our approach per- 
mits us to consider the response of the system to an ac 
field with frequency of the order of several THz which is 
experimentally realizable. 

In conclusion, we have studied the nonlinear response 
of a non-equilibrium quantum dot, subject to a spin bias, 
to an external charge ac field. We computed the spin 
current rectification coefficient of the dot as a function 
of the position of the resonant level and the frequency 
of the applied field. As a function of the resonant level 
position it shows two characteristic sets of three peaks 
(one central and two side peaks) reflecting the spectrum 
of the dot. The heights of the side peaks are consistent 
with the dependence of the relative weight of the reso- 
nant and of the upper level in the density of states on the 
average number of dot electrons. The frequency depen- 
dence of the spin current rectification coefficient shows 
clear evidence of a photon-assisted tunneling - the spin 
is transferred through the dot upon electron absorbing a 
photon to tunnel through the dot barriers. 
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FIG. 1: The spin current rectification coefficient of the dot (in 
units of a 2 e 2 /h 3 ) as a function of the position of the resonant 
level (gate voltage) computed for frequency of the driving field 
Q = O.67 and spin bias voltage eV s = O.57 (solid line) and 
eV s — O.I7 (dashed line). 





FIG. 2: The frequency dependence of the spin current rectifi- 
cation coefficient (in units of /h 3 ) for spin bias eV s = O.I7 
(solid lines), eV s = O.67 (dashed lines) and eV s = 1.2-y (long- 
dashed lines). The resonant level energy is e c = —257 (lower 
panel) and e c = — 57 (upper panel). 




-0.4 1 1 1 1 1 1 1 

10 20 30 

QJy 



FIG. 3: The frequency dependence of spin current rectifica- 
tion coefficient (in units of a 2 e 2 /h 3 ) for spin bias eV s = O.I7 
(dashed lines)and eV s = O.57 (solid lines) and resonant level 
energy e c = — O.27 (thicker lines) and e c + U = — O.27 (thinner 
lines). 



